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Nematic Films and Radially Anisotropic Delaunay Surfaces 
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Abstract. We develop a theory of axisymmetric surfaces minimizing a combination of surface tension and 
nematic elastic energies which may be suitable for describing simple film and bubble shapes. As a function 
of the elastic constant and the applied tension on the bubbles, we find the analogues of the unduloid, 
sphere, and nodoid in addition to other new surfaces. 

PACS. 61.30.-v Liquid Crystals - 82.70.Rr Aerosols and foams - 47.55.D Drops and bubbles 



1 Introduction and Summary 

Composed of elemental air and water with a dash of soap, 
bubbles and foams are fodder for studies in geometry, 
topology, and assembly. Understanding their structure and 
dynamics is essential for many aspects of industrial pro- 
cessing, flow control, and clean drains. The description of 
bubbles follows geometrically simple laws: a single bub- 
ble forms a sphere to minimize its area for fixed volume. 
When multiple bubbles come into contact, their inter- 
faces obey a set of rules, first suggested by Plateau in 
1873 [J and later proven by Taylor in 1974 0. These 
laws, simple and universal, state that at every edge joins 
three faces (at 120°) and every vertex joins four edges 
(at cos _1 (-l/3) « 109.47°). While the angles are set by 
simple force balance, the number of faces-per-edge and 
edges-per- vertex is fixed by area minimization - to change 
these numbers (3 and 4) would require increasing the to- 
tal area. Thus, in some sense, the foam problem is "pure" 
- there are no coupling constants or ratios of parameters 
that determine the ideal bubble shape. In this paper we 
study a simple generalization, in principle realizable in 
soft matter systems: namely, we add in-plane nematic or- 
der to the foam. Are there more general rules for these 
more complicated materials or do the shapes depend crit- 
ically on a new length scale made from the surface tension 
and Frank elastic constants? We conclude that not only 
are the shapes different quantitatively, but that the whole 
"phase diagram" of shapes as a function of applied ten- 
sion is dramatically changed. The nematic order is able to 
support a much richer class of shapes in addition to the 
classical constant mean curvature Delaunay surfaces: the 
sphere, unduloid, and nodoid. 

Ultimately, we would like to predict the scaling be- 
havior of these so-called nematic foams with time as they 
coarsen, in analogy with the standard work beginning with 
Lemlich for normal foams [3,4,5]. There it has been found 
that the average bubble size (R) ~ i 1 / 2 , compatible with 



experimental results. In preliminary experiments on liquid 
crystal foams [6] it was suggested that (R) ~ t 2 . This 
difference from the usual coarsening exponent was sug- 
gested to be a consequence of defects on the surfaces. In 
order to study the nematic foam, we will here concentrate 
on the geometry of general axisymmetric surfaces station- 
ary under a sum of surface tension and nematic elastic 
energies. In a future paper, we will consider the dynamics 
of a double bubble system under diffusion [7] . We begin by 
deriving a functional for this system, then examine single 
bubble shapes, primarily the generalization of axisymmet- 
ric constant mean curvature surfaces, known as Delaunay 
surfaces. 

The functional we consider is simply that of the Frank 
free energy of an axisymmetric area-minimizing surface 
with director field along latitude lines: 



S[r{z)} = 2tt 



dz 



1+ K 



' P 2 

r 

2 



(1) 



We consider the minimizers of this functional as general- 
izations of Delaunay surfaces. The Delaunay surfaces min- 
imize the following functional, which consists of a surface 
tension term and a pressure term: 



S[r(z)] = 2tt 



dz 



1 r 



P 2 

—r 
2 



(2) 



Minimization of @ leads to the Young-Laplace law, p = 
2^H, where H is the mean curvature of the interface. 

Our generalization of Plateau's problem is in a differ- 
ent direction than what has been studied. Unlike the work 
done on surfaces in external fields (e.g. drops and crys- 
tals in a gravitational field) [8"ll9"IITU] , our anisotropy arises 
spontaneously through the nematic director. Further, our 
functional depends on an internal degree of freedom and 
not the extrinsic geometry of the bubble as in the Wulff 
model |llU12j . Previous work on fluid membranes with 
in-plane order have focused on shape transitions [T3TTH] 
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forced in by topology, not curvature. There has been work 
on the Canham-Hclfrich membrane model similar in spirit 
to our work [15] . Other seminal work includes the study 
of nematic order on rigid membranes [16] . Our study of 
bubble shapes complements work on buckling due to de- 
fects [17: and vesicle shape due to tilt order [18] - there 
the energetics and behavior near defects was the focus. 
Here, we find a class of shapes when the curvature energy 
is highly anisotropic. 

In the next section we will expand on our model and 
study an even more general class of free energy functionals. 
We analyze these models and determine whether we have 
nodoids, unduloids, or bubbles as we vary both the elastic 
constant and the applied tension along z. In section 3 we 
describe the actual axisymmetric shapes having reduced 
their profile to tabulated functions. In section 4 we focus 
on the free bubble with no applied tension. We conclude 
with section 5, the conclusion. 



2 The Model 

In this section we put together the requisite energetic con- 
tributions to derive an energy functional for a thin nematic 
liquid crystalline film without any translational order. An 
ordinary soap film is modeled by a surface ScM 3 , which 
minimizes: 

F A = 7 f dA (3) 
Js 

where 7 is the surface tension. With a volume constraint, 
the stationary points of this free energy are constant mean 
curvature (CMC) surfaces which obey the Young-Laplace 
Law, Ap — 2~fH, where Ap is the pressure jump across 
the interface and H is its mean curvature. In the absence 
of any volume constraints we have H — and the surface 
tension does not play a role in the ground state structure. 
As Veit Elser has put it, this problem of minimal surfaces 
is "pure" - it does not depend on any particular interac- 
tions or elastic constants [19]. Here we aim to minimally 
embellish the physics of minimal and CMC surfaces. We 
add a nematic director n, obliged to lie in the surface, and 
include its Frank free energy [20] : 

F= ~ J d 3 x {#1 (V • n) 2 + K 2 [n • (V x n)] 2 

+K 3 [nx (V x n)] 2 } (4) 

where K\, K2, and K 3 are the splay, twist, and bend elas- 
tic constants, respectively. When restricted to a thin film 
of thickness h ~ kT / Ki with planar boundary conditions 
on the inner and outer surfaces, the director will be in- 
dependent of the depth in the film and we may treat it 
as a two-dimensional line field in the tangent plane of the 
surface, S. 

When the nematic director lies in a surface, the twist 
term vanishes; without loss of generality we work in the 



tangent plane to an arbitrary point so that the director is 
n = (n Xi n y ,0), and the normal is in the z direction. Then 

n ■ [V x n] = n x (d z n y - d y 0) + n y (d x - d z n x ) 
+0(d x n y - d y n x ) 
= (5) 

since n is independent of z at this point. Thus the twist 
term vanishes and the bend term in simply (V x n) 2 . 

In addition to the difference between splay and bend 
K\ — K 3 , the surface normal N breaks the bend term 
into two rotationally-invariant pieces. Defining t = n x 
( V x n) , we have two scalars, t n = N • t and t± — |t | — tjf. 
These two types of bend can each have a different elastic 
cost for bend in the surface and out of the surface, K 3 n 
and K 3 ±, respectively. The free energy contribution from 
the nematic is: 

K = \ [ dA {l^ (V • n) 2 + K 3N t 2 N + K 3± tlj (6) 
•* s 

Note that this is, in principle, a functional of both the 
shape of the surface and the director field living on it. We 
will focus on this interplay in the following. 

Recall that Delaunay surfaces are CMC surfaces of rev- 
olution. In order to generalize these for nematic films, we 
will focus on the case where the director field lies paral- 
lel to the lines of "latitude" on a surface of revolution. 
In cylindrical coordinates, (r, <f>, z), these are lines of con- 
stant z. These configurations are automatically splay-free, 
further removing the dependence on K\. In future work, 
we will consider the fusion of two such nematic bubbles 
(i.e. the double-bubble problem [5T]) since this line field is 
consistent with the condition that the director is parallel 
to the Plateau borders (Fig. []} 0. 

In Cartesian coordinates, the director field at (x, y, z) 
is n = (y,—x,0)/r and the normal is N = (x,y,z)/R, 
where r = \j x 2 + y 2 and R = \/ x 2 + y 2 + z 2 . Direct cal- 
culation confirms that V • n = 0. And 

t = n x [V x n] 

Wehave|t| 2 = l/r 2 ,t 2 N = l/i? 2 andi 2 L = (esc 2 8 - l) / R 2 
where 9 is the polar angle in polar coordinates. 

In keeping with the spirit of minimally modifying the 
soap bubble problem, we will set K 3 ± — K 3 n = K 3 
in the proceeding. We shall show that even without the 
anisotropics of K\ and K 3 ± — K 3 n, the shapes of sin- 
gle surfaces are qualitatively modified by K 3 ^ 0. In this 
limit, the nematic free energy only depends on r, the xy- 
radius of the film so that 

Fn = ^ / dA± (8) 

Axisymmetric surfaces may be parameterized by the z co- 
ordinate so that R(z) = [r(z) cos 9, r(z) sin 9, z\. The total 
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free energy as a function of the cylindrical coordinate r(z) 
parameterizing the surface is 



F[r(z)] = 2tt7 J dz ^r\/T 



1 



K 



P 2 

-r 



(9) 



where K = hK^/j with has units of length-squared and 
the where p = AP/j is the reduced pressure. Taking val- 
ues for the lyotropic liquid crystal PAA, K3 ~ 1.2 x 10 -6 
dynes and its surface tension is ~ 35 dynes/cm [22]. Thick- 
nesses of typical soap films are about four microns. There- 
fore, K ~ 1.4 x 10~ n cm 2 . This results in a length scale 
Va ~4x 10~ 6 cm, which is smaller than the thickness of 
the film. We point out, however, that in a dry foam, things 
could be quite different. In particular, since the variation 
of the area is the mean curvature, H, we expect that after 
coarsening, H is close to zero. In this regime, the principal 
curvatures need not vanish and, indeed, our extra term is 
sensitive to them. We will explore this possibility in a se- 
quel 7\ which studies the dynamics of this new class of 
bubbles. In flowing polymer systems or films with mag- 
netic dipoles, K could be considerably larger; in the mean 
time, we focus on this new variational problem and show 
that it leads to novel shapes and geometries, even in this 
simple limit. 

Let's take a step back from our particular model to 
consider a broader class of functionals. We can write our 



energy as: 



F = 2tt7 / dz 



{[l + if 3 (r)]rVl + r 2 -|r 2 } (10) 



We interpret 1 + Kg(r) as a spatially varying surface en- 
ergy. This is unlike the anisotropic surface energies in the 
sense of the Wulff construction for crystal shapes, in which 
the anisotropy refers to a dependence on the direction of 
the normal vector [TT1[T2] . This energy is spatially varying 
in r, and has, to our knowledge, not been considered, in 
contrast to the spatially varying in z case [8ll9lH0] . Note 
that g(r) ought to be even as a function of r, in order 
for there to be the symmetry r —r. For the case of 
latitudinal nematic director lines, g(r) = 1/r 2 . 

We will be comparing the properties of surfaces mini- 
mizing (strictly speaking, extremizing, as we won't study 
the second variation) this free energy with the well-studied 
case of Delaunay surfaces (K = 0). In order to give a more 
complete picture of the properties of surfaces minimizing 
spatially varying surface tension energies, we will begin 
by discussing the axisymmetric minimizers of the general 
family of free energy functionals described by fTO]) . We 
will call these minimizers "radially anisotropic Delaunay 
(RAD) surfaces" . 

Since the integrand is of the form f(r, r), we may take 
advantage of the standard first integrals of the calculus of 
variations. Because / does not depend explicitly on z, (if 
z were time this would lead to conservation of energy) , we 
have: 

f~r% = T (11) 



di 



where T is a constant. This last equation may be inter- 
preted as a conservation of force along the z-axis, with T 



an applied tension or compression, depending on its sign. 
For a free bubble, T = 0, which is the analog of the sphere. 
Setting this to other values will yield RAD surfaces anal- 
ogous to the other Delaunay surfaces, such as the nodoid 
and unduloid. We shall see that the types of RAD surfaces 
that are possible as we vary T will depend on the value of 
K in a nontrivial way. 

If we insert our expression for f(r,r), we can rewrite 
CD: 



±7 



2nr 



VI+7 2 



[l + Kg(r)]=T + nr 2 AP (12) 



This is an equation for force balance in the z-direction 
0: on the right hand side, the two terms correspond re- 
spectively to the applied tension and the pressure-induced 
force on a surface with xy cross-section 7rr 2 . Likewise, the 
left side is the effective surface tension 7[1 + Kg(r)] mul- 
tiplied by the length of the boundary, 27rr. The factor of 
±(1 + f 2 ) -1 / 2 projects the z-component of this force. The 
ambiguity in sign arises from the fact if r passes through 
infinity, the sign of the surface tension force must change 
as well. In subsequent analysis we will need to be careful 
about which branch of the radical we choose, in particular 
in the case of the nodoid where r does indeed goes through 
infinity and the z-component of the tangent vector to the 
curve changes sign. We will cross that cut when we come 
to it. 

By rewriting for AP, we can derive a generalization of 
the Young-Laplace rule: 



AP 



2 7 [1 + Kg(r)} 
ryl + f 2 



T 

(13) 



Choosing the film to be symmetric about z — 0, we know 
that r(0) = and if we define R = r(0), the radius of the 
"waist" then 



(14) 



When K and T both vanish, this reduces to the standard 
Young-Laplace rule. 

En passant we note that in the case of contact between 
two bubbles, the effective surface tension 7[1 + Kg(r)] is 
identical for the three partial caps shown in Fig. [1] since 
the value of r is the same where they meet. As a result, the 
triple junction will join the three surfaces at 120°, just as 
in regular foams. However, the equality of surface tensions 
is an accident of axisymmetry and in general, Plateau's 
rules will not hold. 



2.1 Families of RAD surfaces 

In this section we analyze the families of RAD surfaces as 
we change T and K. First, let us measure T in units of 7 



We note that in the study of CMC surfaces, this is an ex- 
pression for the magnitude of the weight vector |23| , which can 
be used to show that CMC surfaces must asymptote to Delau- 
nay surfaces at infinity. Whether this notion can be generalized 
to our problem is an open question. 
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K = 0.01 




Fig. 1. Double RAD bubble, with g(r) = 1/r 2 and K = 0.01. 
Three T — caps meet at 120°. 



as we did for K and p, i.e. T = T/(yvy). The force balance 
equation with our radial anisotropy function specialized 
to g(r) — 1/r 2 (the nematic-derived case) now reads 



l + K 



P 2 

r 



T 



(15) 



If P 7^ C0; w e m ay also scale out the dependence on p via 
f = pr/2, z = pz/2, f = pT/2, and K = p 2 K/A: 



K 



f 2 =T 



+ r 



r2 



(16) 



In the following analyses, we will drop the tilde from our 
variables unless there is any ambiguity. There remains 
only one free parameter, K. 

We begin the study of the family of RAD surfaces with 
a given K by viewing T = T(r, r). In Fig. [2] we have plot- 
ted the embedded surfaces [r, f, T(r, f)] for both branches 
of the radical. Contours of constant T correspond to val- 
ues of r, r on different axisymmetric surfaces with a fixed 
tension. 

We start by considering the positive branch phase por- 
trait first - the features sketched below should be consid- 
ered as a taste of things to come in our later analysis. The 
diagram of this branch of the T-surface represents points 
on the generating curves of RAD surfaces where the tan- 
gent vector points upwards (has a positive z-component). 
In Fig. [3] we have set K = 0.01 and find that there are 
several families of contours which cut across this branch 



2 We have not examined the p — surfaces, but these will 
be to minimal surfaces as general RAD surfaces are to CMC 
surfaces, and will hence be a generalization of catenoids. 




Fig. 2. The surface T(r,f), when K = 0.01. Each connected 
constant height contour (green) corresponds to r, r for a gen- 
erating curve for a RAD surface. 

T(r,t); K=0.01 



t o 




Fig. 3. The constant T contours on the positive branch of 
T(r,r). Closed contours correspond to unduloids, the single 
maximum is a cylindrical solution. 



of the T-surface. First, the T = contour represents the 
curve (r, r) for a free bubble, indeed, there is a range of 
T such that the contours imply that r — ±oo as r — > - 
these all represent closed bubbles. Note that among these 
there are those which remain completely on the positive 
branch (convex bubbles) , and some which run to the neg- 
ative branch (Fig. [4]) as well (concave bubbles). There are 
also closed curves with T > on the positive branch which 
represent solutions periodic in z. These are the so-called 
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T(r,t);K=0.0I 




o.o 



0.2 



0.4 



0.6 



0.8 



1.0 



Fig. 4. The constant T contours on the negative branch of 
T(r,r). Each contour joins with one on the positive branch as 
it passes through r = oo. 



unduloidal solutions. The single point solution is a cylin- 
der (since r =const). Finally, there are also curves with 
T < which "jump" to the negative branch and back, 
that is r — ±oo for nonzero r - these are the nodoids, 
the periodic self- intersecting surfaces, where the tangent 
vector turns in a loop in each repeat unit. Note that not 
all surfaces cut the negative branch of the T-surface, but 
all must cut the positive branch. To see this, note that 
in (fro]) when T > 0, we must choose the positive branch, 
while when T < we have solutions for both signs of the 
radical, in particular for large r. 

This can be generalized for any g(r). Contours of con- 
stant T are either closed (unduloids), run off to r = ±oo at 
r — (bubbles), or run off to f = ±oo at finite r (nodoids). 
Here we sketch how to partition the full space of RAD 
surfaces into these three classes when g(r) is arbitrary by 
studying the special case g(r) = l/r 2n for n G Z + . This 
corresponds to assigning an energy to the lines of latitude 
of K 2n , where k is their curvature in R 3 , and one of the 
principal curvatures of the surface. We label the second 
argument of T(-, •) as t = r in the following and will study 
the function T(r,t). 

It's a fact from the theory of bifurcations in dynami- 
cal systems that the critical points of the T(r, i)-surface 
indicate where the geometry of the contours change, and 
hence those values of T for which the geometry of the 
RAD surfaces change. For example, near a local minimum 
or maximum, the contours of constant T are closed curves 
encircling the extrema. However, at a saddle point of the 
T-surface, these closed contours are separated into two 
groups; the separatrix is a self-intersecting curve and cor- 
responds to a singular topological change in our surface. 



Thus, if we find the critical points and their types (min- 
ima, maxima, saddle) , then each of the families mentioned 
will correspond to families of contours which lie on points 
of the surface between the critical points. In order to find 
the critical points, we calculate the gradient of T and look 
for its zeros. 



VT = ± 



1 - K{2n - 1)7 



VT+t 2 



- 2r,- 



rt(l + Kr~ 2 
(l + t 2 ) 3 / 2 



When VT = 0. 



± 



1-K(2n- l)r~ 2n 

a/T+1 2 

rt{l + Kr- 2n 



= 2r 
= 



(17) 
(18) 



Since r is real and positive, we must have t = and so all 
critical points of T must lie on the line t = 0. The Hessian 
H provides information on whether a critical point is a 
saddle, minimum, or maximum. At t = 



H 



±2Kn(2n 



- l) r - 2n - 






=pr(l + Kr 



-In- 



(19) 



Since H is diagonal, it is easy to see that the directions 
of maximum change in T are along the r and t directions 
with the corresponding eigenvalues on the diagonal of the 
matrix. Note that on cither branch, H t t does not change 
sign. Thus on the positive branch it is always negative and 
the t = line is ridge- like; similarly on the negative branch 
the t — line is valley-like. Thus our work is simplified 
since we need only check H rr to determine which type of 
extrema we have. 

As we have discussed, the T = solution has signif- 
icance as the "free" solution with no imposed tension at 
the boundaries, but is also important in another respect. 
Solutions with T > cannot have portions on the nega- 
tive branch, as then both terms on the left hand side of 
the force balance equation (fT3|) would be negative. This 
implies that the T = solution lies at the boundary be- 
tween surfaces which may be parametrized simply as r(z), 
and those which cannot - a result of the behavior of the 
z-component of the tangent vector. Luckily, a half-repeat 
of any surface can be parametrized as z(r), which saves 
us from having to look at a completely parametric repre- 
sentation (r(s),z(s)). 

Along the ridges/valleys we have 



T = ±r 1 



K 



or 



(r 2n + Kf - r 



4n-2 



r 2 ) =0 



(20) 



(21) 



We denote r to be the value of r(z) when dr/dz = t = 0. 
All possible values of rg are real roots of (|2ip. Letting 



u = r , we see that the real roots rj of the degree An + 2 
polynomial in r arise from positive real roots Uj of the 
degree 2n + 1 polynomial in u 



P{u) = (u n + K) 2 



,2n- 



1( T . 



(22) 
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On the T-surface, the points (rj, 0) are where the contour 
corresponding to a RAD surface intersects the r-axis. Note 
that when T(rj, 0) is a critical point, rj must be a double 
root of P(r 2 ) = 0. We will use this fact many times in 
what follows. Though it would be tempting to associate 
each distinct positive root Uj for fixed T with a distinct 
RAD surface, a single generating curve may have multiple 
points where dr/dz = 0. 



2.2 Reduction to Quadrature 

Indeed, while it's simple to examine the picture of the T- 
surface to see which points are part of the same contour, 
the mathematics is more subtle. To flesh this out, we take 
a small detour to integrate our force balance equation by 
quadratures for the function z(r). Again, there is a sign 
ambiguity in the integrand, the same sign ambiguity that 
we have seen throughout. We now set the sign by the 
boundary condition that ^-component of the curve's tan- 
gent vector is positive where we begin the integration. We 
choose the origin along the z-axis at the place where f = 0. 
At this place P(u) vanishes, by construction, and we have 



z(r) - 



(T + u) u n - 1 du 



r2 W 2^J (u n + K) 2 - u 2 "- 1 (T + uj 



= (23) 



Note that the denominator vanishes as yfu near the lower 
limit and thus the singularity is integrable. Also, since 
we only integrate between consecutive roots of P(u), the 
sign of the integrand remains fixed. Of course there are 
several roots of P(u). Let Uj for 1 < j < 2k + 1 be a 
positive real root of P(u), with = uo < U\ < ■■■ < 
U2k+i-, so that 2k +1 is the total number of real roots (an 
odd degree equation with real coefficients is guaranteed 
an odd number of real roots since complex roots must 
come in conjugate pairs). The values of T, K control the 
coefficients of P(u) and hence whether k = 0, 1, . . . , n. The 
integral for z(r) can only remain real when P(u) > and 
the intervals of u such that this is true are also bounded 
by the Uj. Therefore, only for r 2 such that P(r 2 ) > will 
the integral from Uj to r 2 result in a real valued z(r). 

Each interval over which P(u) > corresponds to a 
distinct type of RAD surface. When K = 0, the only 
closed bubble solution is the sphere when T = 0. How- 
ever, when K ^ 0, P(0) = K 2 > for any T and the first 
valid interval is r 2 S (0, Ui). This corresponds to closed 
bubbles, as r — means that the generating curve touches 
the z-axis. The existence of this interval for K there- 
fore means that there will be a closed bubble solution for 
any T. 

To exercise our formalism, we review the well-studied 
K = Delaunay surfaces. Taking K — in our equations 
directly yields 



P(u) = (u n + 0) 2 - u 2 ™- 1 (T + u) 



,2n-l 



u~{t + u y 



(24) 



We will drop the factor of u 2n ^ 1 as it yields 2n—l degen- 
erate RAD "surfaces" u — 0. Then the system is indepen- 
dent of n, and we recover the classical Delaunay surfaces. 

In Fig. [5]we plot T(r, 0) versus logr; this is a plot 
of the cross section of the T-surface through the plane 
t = 0, flipped on its side. This graph is a "tomographic 
slice" of the a general T-surface like that in Fig. [2] except 
with K = 0. The two branches of the boundary curve 



are T(r, 0) = ±r — r 2 or, in other words, u±(T) 



: (T) 



are the solutions to the quadratic equation P(u) = u — 
(T + u) 2 = 0. The region between the two curves satisfies 
P(u) > 0, and has been shaded blue and green, for the pe- 
riodic surfaces and the closed bubble surface, respectively. 
The generating curves for Delaunay surfaces are found by 
setting K — in (|2lij) . so they are described by 



z(r) 



(T + u)du 



2^fu^Ju — (T + u) 2 



(25) 



where u_ is the smaller root of P(u) 



u± = ]-{l-2T± Vl -AT 



(26) 



When T < j, the roots are real and u± > 0. Moreover, it- 
vanishes only for T = and so the only closed surface at 
K = is the free, tensionless bubble. One can see easily 
from z(r) that it is a sphere in this case. Recall that this 
solution with T = separates surfaces with contours that 
lie entirely on the positive branch of T from those that do 
not. For the Delaunay surfaces, the difference between the 
unduloids (T > 0) and the nodoids (T < 0) is straightfor- 
ward: the nodoids all self-intersect, and the unduloids do 
not as shown in Figs. [6] and [3 

When the two roots coincide, the integrand diverges 
over the integration region unless r 2 is constant and equal 
to u±, i.e. a cylinder. The discriminant of the quadratic 
equation, Q(T) = 1 — 4T vanishes at this point, and one 
can also check that T = 1/4 is the only critical point on 
the T(r, t) surface. Since points with t = on the positive 
branch of T are on a ridge, the cylindrical solution is the 
maximum of T. 

Finally, we may study the "period" of the Delaunay 
surfaces as we change T as a further diagnostic of the 
resulting bubbles. We define the period as the twice the 
height between the endpoints of a surface or Z = 2[z(r + ) — 



z(r_)], where r 2 = 



u + (T + u)du 

_ ^/u~P(u) 

u + (T + u)du 

_ ^u(u — u-)(u + — u) 



2T 



-K 



2r,E 



U- 
14-1 



(27) 
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™y\ : 




jr cylinder" 


■ ■ ■ 


J free bubble 



-0.3 -0J2 -0.1 



0.0 
T 



o.i 



0.2 0.3 



Fig. 5. A log-plot of the t = slice through the T-surface for 
K — 0. The shaded portion of the graph maps out the range 
of u — r 2 (z) for each Delaunay surface as a function of T. 



valid when T < 1/4 and where K(-) and E(-) are the com- 
plete elliptic integrals of the first and second kind 0. As 
shown in Fig. Q3J as T — ► — oo, the period goes to zero 
(corresponding to tightly wound nodoids), and at T = 0, 
U-/u+ = 0, r + = 1 and the period, and hence the diam- 
eter of the sphere is 2. As T — ► 1/4, the period goes to 
7r. This might be unexpected, as the period of a cylinder 
ought to be undefined, or perhaps zero or infinite. How- 
ever, this result just implies that the cylindrical limit is 
reached by the vanishing of wiggles in the surface, rather 
than singular behavior in the period function. 

Having cut our teeth on the K = case, we move on 
but refer the interested reader to other sources 24.25 . As 
we mentioned, an essential difference between the K = 
and K ^ case is the possibility of additional distinct 
solutions for fixed T. If K ^ 0, we can find new solu- 
tions as follows: since P(u) > in the interval (0, U\), and 
P(ui) = 0, it follows that P(u) < for u £ (u l7 u 2 ). Thus 
for r 2 £ (u 2 , 1*3), P(r 2 ) > and we find a new valid solu- 



dO 



o \/l - k 2 sin 2 9 
sin0 dz 



F(<p,k) = 



i/(l- z 2 )(l- k 2 z 2 ) 

E(<p, k) = f Vl-k 2 sin 2 0d6 
Jo 

f 8 ^ VI - k 2 z 2 dz 




Increasing T 

Fig. 6. Generating curves of Delaunay surfaces. Nodoids are 
red, unduloids are green, and the spherical solution (T = 0) is 
in black. Dashed curves show a few repeats of a nodoid and an 
unduloid. 



K = 0;T = -1.0 




and K{k) = F(ir/2, k) and E(k) = E(tt/2, k) 



Fig. 7. Comparison of a Delaunay unduloid (a) and nodoid 
(b), at T = 0.2 and T = — 1 respectively. Note that the loops 
formed by the self-intersections in (b) point "inwards" , so that 
this is an example of an "innie" nodoid (c.f. Fig. I17H - ). 



tion of the force balance equation. We may evaluate the 
integral for z(r) either from r 2 to U2j+i or from U2j to r 2 - 
the difference will be whether we define the plane z = at 
the "outer radius" or at the "inner radius" of the surface 
and corresponds only to a shift along z. Continuing with 
this reasoning we will find distinct surfaces by integrating 
values of r 2 £ (u2j, u 2 j+i) for j = . . . , fc, so that we have 
k + 1 distinct surfaces. 

As mentioned earlier, there is a connection between 
multiple roots of P(u) and critical points on the T-surface. 
If there are repeated roots, then by perturbing T slightly, 
we can interpret the solutions as lying on a separatrix. 
Depending on the parity of the colliding roots, either two 
intervals of positive P(r 2 ) join (when u 2 j — > u 2 j+i for 
some j) and two families of surfaces join into one, or an 
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interval where P(r 2 ) > vanishes and the limits of the 
integral degenerate into a point and we have an infinite 
cylinder of radius r 2 = U2j+i = "2j+2 for some j. 

Therefore, studying the discriminant of P(u), which 
we call Q(K,T) will reveal how many positive real roots 
P has, and hence how many distinct RAD surface^. For 
simplicity, we specialize to n = 1, where 

P(u) = (u + K) 2 - {T + ufu 
Q(K, T) = -(K- T) 2 (4K + 27K 2 - 18KT - T 2 + 4T 3 ) 
= -{K -T) 2 Q 3 (K,T) (28) 

where the last equation defines the cubic factor Q 3 (K, T) . 
Since P(u) is cubic, it may have one or three real roots. 
When the discriminant vanishes, roots collide. This may 
happen if K = T or at a root of Q 3 . K — T implies 
that the roots of P(u) are {1, -K, ~K}. When K > the 
repeated root corresponds to non-physical, imaginary val- 
ues of r and so the vanishing of Q(T, T) is of no particular 
physical consequence there. We will discuss the K < 
case only briefly at the end of this section. We thus focus 
on the roots of Q 3 (K,T) which are precisely the criti- 
cal points in T, the values of tension at which there is a 
crossover between a system where there is one bubble (one 
real root), and one where we have a bubble and another 
surface (three real roots). We'll discuss the geometry of 
these surfaces later. 

Moreover, the number and nature of the critical points 
of T varies as we changed K. Note that as K changes, the 
positions of the roots of Q 3 (K, T) move as well. Hence we 
may expect that if we tune K through certain critical val- 
ues, the critical values of T will themselves collide or split, 
and so we gain or lose entire families of RAD surfaces. It 
will therefore be of use to consider the discriminant of 
Q 3 (K,T). K — will necessarily be a critical value in K 
because there the stationary surfaces are simply the or- 
dinary Delaunay surfaces along with multiple degenerate 
u = roots in P(u) which we dropped in our discussion 
above. Since this case is understood, we focus on nonvan- 
ishing critical values of K . 

Since the relevant factor of the discriminant Q(K,T) 
is cubic, we may have one or three (real) critical points in 
T. The discriminant of Q 3 (K,T) is 

D 3 = -16K(-1 + 27K) 3 (29) 

Descartes rule of signs tells us that Q 3 (K,T) will always 
have one (real) negative root in T, which we label T\ . If the 
other two roots are complex, then let T 2 — T 3 * = /i + iv. 
Since D 3 = 64{T 1 -T 2 ) 2 {T 1 -T 3 ) 2 {T 2 -T 3 ) 2 we see that in 

the case of one real root D 3 = —256 [(Ti — ^) 2 + v 2 ] 2 v 2 < 
0. When all three roots are real, however, D 3 is real and 
positive. Thus when < K < 1/27 there are three critical 
points, but for larger K there is only one critical point on 
the T-surface. 

4 Recall that for a general cubic, ax 3 + bx 2 + cx + d with roots 
ri, r-2, and rz, we have rxr^ri, = —d/a, rir2 + r2r3 + rirz = c/a, 
and n + r 2 + r 3 = -b/a. It follows that b 2 c 2 - 4ac 3 - 4& 3 d - 
27a 2 d 2 + I8abcd = a i (r 1 - r 2 ) 2 (ri - r 3 ) 2 (r 2 - r 3 ) 2 . This is the 
discriminant of the cubic. 



K=0.137 




T 



Fig. 8. The t = slice through the T-surface at K = 0.137. 
Only two regimes are observed here; nodoids and bubbles for 
T < Ti and bubbles only for T > Ti. 



K=l/27 




T 



Fig. 9. The t = slice through the T-surface at the critical 
point of K = 1/27. Compare Figs.|8](A" > 1/27) andflOlfA' < 
1/27). 

Having deduced so much with so little, let's unravel 
what this implies for P{u) in each case. We begin with 
the case K > 1/27, as T\ < is our only critical point 
and the situation is somewhat simpler. When T < T\, 
there are three real roots and we will have two intervals of 
positivity, (0, u\) and {u%, 113), and hence we have two dis- 
tinct RAD surfaces, a closed bubble, which we call a "tiny 
bubble" and another surface we call a nodoid in analogy 
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K= 10" 



Q 



-10 



-15 



unduloids - 




"innie" nodoid + 
tiny bubble 



-0.4 



-0. 



0.0 

T 



0.2 



0.4 



Fig. 10. The t = slice through the T-surface with K = 10~ 4 
(qualitatively similar to Fig. [2] where K — 10 -2 ). There are 
four regimes, nodoids and bubbles for T < Ti, bubbles only 
for I\ < T < T2, unduloids and bubbles for T 2 <T <T 3 , and 
bubbles only for T > T 3 . 



with the Delaunay surface since T < 0. T\ corresponds 
to the point where the nodoids join with the tiny bubble 
(i.e. u\ collides with u 2 and they move off the real axis). 
When T > T\ there is one real root of u and we have one 
interval for which P{u) > 0, and the only valid integra- 
tion interval is the one corresponding to the closed bubble, 
namely (0, ui). We note that this root is a continuation of 
u 3 from when T < T x . In Fig. K = 1/27 + 0.1, we show 
a semilog plot of the t = slice of the T-surface depicting 
the space of surfaces. 

On the other hand, if there are three critical points 
on the T-surface (when < K < 1/27), things become 
more complicated. Let the three critical values be then 
T x < < T 2 < T 3 . When T < T x , P(u) has three real 
roots, and so we have two distinct RAD surfaces, a "tiny 
bubble" and a nodoid again. T\ corresponds to the point 
where the nodoids join with the tiny bubble. For T\ < T < 
T 2 , we have only a closed bubble solution, as now P(u) 
has only one real root. We'll sometimes call this the large 
bubble, to distinguish it from the tiny bubble solutions, 
though we will see (Fig. [TTj) that they may be taken to 
each other via a non-singular path surrounding the point 
(T = 1/3, K = 1/27). 

We have seen that for T < 0, there is always a portion 
of the contour on the negative branch of the T-surface, so 
the large bubble with T < must be concave at r = 0, 
which we can verify by calculating d 2 z/dr 2 at r = 0. For 
T > 0, the bubble is convex, and T = is the free bubble 
case which we will consider in much more detail in section 
4. In that case, we must consider d A z/dr A at r = to 



self- intersecting 
closed bubble 



non-self-intersecting 
closed bubble 




•1.0 -0.5 

/ 

closed bubble 



cylinder + 
tiny bubble 



tiny bubble 



establish convexity. 



Fig. 11. RAD surface phase diagram, when g(r) — 1/r 2 . The 
solid red curve and dotted dark blue curves separate regions 
with one RAD surface from regions with two RAD surfaces 
and are the curves defined by the vanishing of Qs(K,T) (eqn. 
I28[) . The period / height function for the surfaces diverges to 
+00 for T < and —00 for T > on the solid red line. The 
dotted green line is the line T = K. Together these curves are 
where the discriminant Q(K, T) = 0. The solid blue and dotted 
orange lines are calculated from the vanishing of the height 
Z X {T,K) and period Z 3 (T,K) feqns. l33l32ll , respectively. Note 
that the non-self-intersecting closed bubble solutions can be 
divided into bubbles convex at r = (T > 0) and concave at 
r = (T < 0). 



For T2 < T < T3, P(u) has three real roots again, and 
we have both tiny bubbles and periodic non-self-intersecting 
RAD surfaces analogous to unduloids. At T3, the undu- 
loids degenerate to a cylinder, though the tiny bubble is 
still there, and past T3 there is only a tiny bubble solu- 
tion. In Fig. [10] we show a semilog plot of the t = slice of 
the T-surface for this case. Fig. [9] shows the limiting case 
K = 1/27, which behaves most like the K > 1/27 case. 
Since T 2 = T 3 , the unduloids in the < K < 1/27 case 
are not visible. Therefore the change in character of the 
RAD surfaces as a function of T is controlled only by the 
critical point at T\. 

In Fig. [TTJ we show a "phase diagram" in T and K 
of the types of surfaces possible. From this, we can see 
the "critical points" at (T, A") = (1/3, 1/27) and (0, 0) are 
where certain "first order" (branch cut) lines end. We con- 
clude this section with a few sentences on the unphysical 
but mathematically enlightening case of K < 0. When 
K < 0, Q3 only has one real root Ti (actually the con- 
tinuation of the root T3 when K > 0), which lies between 
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the regime of unduloids and tiny bubbles and the regime 
of closed bubbles. However, now the (K — T) factor of 
the complete discriminant Q of P(u) becomes relevant. It 
defines a curious boundary between the nodoid and the 
unduloid. Though the line T = K separates two regions 
for which there are two RAD surface solutions, the tiny 
bubble and a periodic surface (nodoid when T < K, un- 
duloid on the T > K side), on the line itself, the tiny 
bubble caps off the periodic surface to form only one so- 
lution, a larger closed bubble. This behavior generalizes 
the nodoid to sphere to unduloid transition at T = for 
the Delaunay surfaces. This observation thus unifies the 
periodic surfaces somewhat. Below we will again concen- 
trate on surfaces with K > 0, though our formulas are 
applicable to all K. 



3 Shapes 

We finally turn to the particulars of the allowed shapes as 
we explore the contours of constant T. When n = 1 the 
radicand in (J22J) is cubic in u and thus the integral may 
be evaluated in terms of elliptic functions: 



z 3 (r) 



(T + u)du 



(30) 



T + m 
%/w 3 = u 



:F 



U 3 



V>3 — U 2 



u 3 - u 2 

u 3 - Ul 



+Vu 3 — u\E 



U 3 



U3 - U 2 



u 3 - u 2 

U 3 - Ul 



where -F(</>, k) and E(<j>, k) are the elliptic integrals of the 
first and second kind, respectively, and the Ui are the roots 
of P. The domain of this function is u 2 < r 2 < u 3 when we 
are looking at nodoids and unduloids. For the large bubble 
case, P(u) has only one real root u\, so we must break this 
naming convention when interpreting the formula above 
by taking Ui, u 2 to be the complex conjugate roots and u 3 
to be the real root. In this case, the domain is < r 2 <u\. 
The tiny bubble shape can be expressed in this form as 
well: 



zi(r) 



(T + u)du 
~2^JpW 



(31) 



T + Ul 

V u 3 - Ul 



F 



>u 3 



"2 



Ul 



+ V U 3 — UiE 



iu 3 



u 2 — r" 



u 3 - u 2 
u 3 - Ul 



(u 3 — r 2 )(ui — r 2 ) 
u 2 — r 2 



Here < r 2 < u\. The integral could describe the closed 
bubble as well, but our choices of elliptic integrals here 



turn out to be on the wrong branch for that case. For 
large \T\, the root of P(u) is u x - T ^ 2K ~ K 2 /T 2 . Since 
in this limit u± is small, the remainder of the radicand 
varies slowly in u € [0, Hence in this limit the bubble 

T I 

pf| V u l 



, the tiny 



shape becomes spherical, Zi(r) 
and closed bubbles become spherical as T — * oo. 

In Fig. [12] we show some representative generating 
curves in the case < K < 1/27, specifically with K = 
0.01. Note that as T increases we go from branches of 
nodoids which begin at z — and end at some z ^ 0, to 
large bubbles which end at r = 0, and finally to branches 
of unduloids which also end at z ^ 0. The generating 
curves for K = 0.55 > 1/27 are depicted in Fig. [T3l where 
we transition from nodoids to concave to convex bubbles 
as T grows. We may also analyze the behavior of the solu- 
tions as a function of T by considering the period function 
as we did for Delaunay surfaces: 



Z 3 = 



(T + u)du 



T 



yj{u- Ui){u - U 2 )(u 3 - u) 
+ Ul 



(32) 



\fu 3 - u 



--K 



u 3 - Ul 



u 3 — u 2 



-2\/u3 — uiE 



u 3 - Ul 
u 3 - u 2 



Z 3 only yields the periods for nodoids and unduloids, and 
the complete elliptic integrals are chosen so that their 
moduli satisfy < k < 1 when P(u) has three real roots, 
so that the integral is real. For the tiny bubble, we have: 



Zi 



(T + u)du 

7W) 



2zi(0) 



(33) 



In the case of three roots Zi also gives the height of the 
closed bubble. In the large \T\ limit, since these bubbles 
become spherical, Z x ~ 2ri = l^fui ~ —2T/K, with the 
sign chosen to agree with that of Zi . The notations for Zj 
and Zk are independent of each other, with j referring to 
the number of real roots required to exist for Zj to make 
sense, and k referring to the upper limit of integration Uk ■ 

In Fig. [M] we show the behavior of Z 3 and Zi for var- 
ious values of K . The K > period functions Z 3 and Zi 
both have divergences at T = T\ and T — T 2 ; this is visible 
in the plot as vertical lines where three curves Z 3 (T < Tj), 
Zi(T <Tj), and Z%(T > Tj) join. On the phase diagram 
(Fig. [TTj) . the divergences lie on the solid red line; the di- 
vergence to +00 at 7\ < and — oo at T 2 > 0. Physically, 
this signals a loss of mechanical stability in these systems, 
as an "elastic modulus" of dF z /dz = dT/dZ goes to zero. 
From our discussion earlier, T = Ti is where the nodoid 
and tiny bubble solutions merge. Fig. Qj)] illustrates this 
for K = 0.01 and K = 0.55. We can also think of this 
process as the tiny bubble "tearing off" as we increase the 
compression on a large bubble (decrease T) past Ti . 

The nodoids close to Ti have the "loops outside" as 
shown in Fig. fTTH . as opposed to the usual "loop inside" 
nodoids in the Delaunay surface case shown in Fig. [7Jd. In 
terms of Z 3 , we see that "innies" are when Z 3 < and 
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"outies" are when Z3 > 0. The period diverges at T = T\ 
in Z 3 , which appears for all K > but not K = 0. The 
fact that the divergence is positive guarantees a value of 
Ttoms < Ti < where Z3 = and a torus appears. We 
show the tori at K = 0.01 and K = 0.55 in Figs. E^f. 
The curve T torus (K) is also drawn on the phase diagram 
(Fig. [n} as well. Note that T torus for If = 0.01 is very 
close to T\ and the generating curve for the K — 0.55 
torus is quite close to circular. In fact, we find numerically 
that as if — > 00, Ttorus ~ —3K — 2/3 and the generating 



K = 0.0I 



curve approaches ± 



a circular torus 



with ring radius y3K and tube radius 2/3. These length 
scales can be derived analytically by finding the first few 
terms in a series solution (in \[K) for 112 and U3, assuming 
the empirically derived behavior ofT iorus . 

Finally, we turn to the closed bubble solutions: as T 
grows from negative to positive, the bubbles change from 
concave to convex at r = 0, this follows from the expan- 
sion z\{r) — jj^(ui~ r 2 ) + 0(r 4 ). From the divergence in 
Z\ at T\ and T2, we know that the height of the bubble 
goes to infinity when it joins with the nodoid and undu- 
loid. Closed bubbles with Z\ = are tangent to them- 
selves at the origin, and so Z\ (T, K) = forms the line 
between self-intersecting and non-self-intersecting bubbles 
(Fig. fTTj) . We show concave bubbles with K = 0.55 in Figs. 
[TTb.c, and a bubble close to T 2 with K = 0.01 in Fig. [PTa. 
The splitting of such a protruding bubble into a tiny bub- 
ble and an unduloid is depicted in Fig. [TBI when K = 0.01. 
We add that this shape is reminiscent of those found in 
Canham-Helfrich membrane models for vesicles [26] . 

Note that although the unduloids end at T = T3, there 
is no indication in the period function; indeed, as with 
the pure Delaunay case, the cylindrical solution actually 
has a finite period as a limit of unduloids. Neither of the 
elliptic integrals making up the period function has any 
divergence at this point. Hence, in contrast to the behavior 
at T\ and T2, the tiny bubble solution at T3 is completely 
regular. 



4 Free Bubbles 

From this point forward we discuss only the T = case 
of free bubbles. In future work we will start from these 
equilibrium shapes to study the dynamics of the nematic 
double bubble [7]. As before, we define z = at the "equa- 
tor", where dr/dz = 0. We define R = r(0) as the waist 
of the bubble. Specialized to T = and g(r) = 1/r 2 , ([14]) 
reads 



P = 



R 



1 



R? J 



(34) 



As we argued, when T = there is a closed bubble for 
all K. Accordingly, R will be analogous to the radius of 
this bubble. We define here the new dimensionless quanti- 
ties K = -£5 and p = pR = 2(1 + K). Our previous choice 
for dimensionless quantities had the advantage of remov- 
ing p, but had the disadvantage that the K — > 00 limit is 




Fig. 12. Generating curves (for the nodoids, large bubbles, and 
unduloids) at K = 0.01. The surfaces are color-coded: nodoids 
are red, concave closed bubbles are purple, T = is black, 
convex (at r = 0) bubbles are blue, and unduloids are green. 
Dashed curves show a few repeats of an "innie", an "outie", 
and an unduloid. 



obscured; this turns out to be an important limit for dy- 
namics. Our choice here represents a change in viewpoint. 
Instead of viewing p as a given parameter that we can re- 
move via scaling, p is now a quantity which we calculate 
given R and K . These correspond to constant pressure ver- 
sus constant size approaches, respectively. In these units, 
the first integral (flU)) determines the bubble shape via: 



2 z 
~R 



pudu 



a* y 4(u + K) 2 -~p 2 u 3 
1 (1 + K)udu 

{u + k) 2 - (1 + k) 2 u 3 



(35) 



We can factor the cubic in the radical to find that the 
only real root is at u = 1. This is consistent with our 
more general arguments (Fig. Ill[) : when there are three 
critical points in T, T — lies in the large bubble region, 
while when there is one critical point T\ < 0, = T > T\ 
is in the closed bubble region. Written in terms of the 
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Fig. 14. The functions Zz and Z\ , which give the period of unduloids and nodoids and the height of closed bubbles, respectively, 
(a) The period surface Z(T,K). The Z3 surface is red and the Z\ surface with K < 1/27 is blue, and Z3 with K > 1/27 is 
green. The single period function Z at K — is the dashed and purple edge of the red Z% surface. The phase diagram (Fig. Hip 
gives a "top view" of this surface, (b) Constant K slices through the period surface. Again, Z3 curves are red and orange, Z\ 
curves are blue for K < 1/27 and green for K > 1/27, and K = is the dashed purple curve. 



K = 0.01, K = 0.55,T = -1.2 K = 0.55, T = -1.85 K = 0.01, T = -0.209 




Fig. 17. From left to right: (a) A closed, protruding bubble at K = 0.01 with T 2 - 10~ r , where T 2 = 0.1887. Concave bubbles 
at K = 0.55 with T = -1.2 (b) and T = -1.8 (c) . (d) An "outie" nodoid at K = 0.55 with T = -1.85; Ti = -1.848 at this 
K. Tori at K = 0.01 (e) and K = 0.55 (f). For K = 0.01, T torus = Ti - 4.48 x 10~ 8 where Ti w -0.209 and for K = 0.55, 
Ttorus = -2.232 and Ti = -1.848. 



elliptic functions F and E, we have: In Fig. [T5] we show the generating curve for several val- 

ues of K. Note that the shape of the curve becomes more 
z ( r ) _ u - p I / \ ~ W / \ ~ u + I oblate as K increases. Intuitively, this makes sense, as the 

R y/1 — u- y y 1 — u+ y 1 — U- J K/r 2 term in the energy penalizes surfaces with more ma- 

terial close to the axis. 





/'- 


u + \ 








1- 




1-U+ 




U+ ' 


Vi-«- 



+^nE \Jt-S.JZ-±t (36) 



u± = ^—(2 + k±i\/4:K + 3K 2 

2(1 + if) 2 V 
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Z 




T — is black, convex (at r = 0) bubbles are blue. Dashed 

curves show a few repeats of an "innie" and an "outie". Fi g- 16 - The transition between large bubble to unduloid and 



tiny bubble for K — 0.01. Note the divergence again in the 




(c) (d) 

Fig. 15. The transition between nodoid and tiny bubble to 
large bubble for K = 0.01 ((a) and (b)) and K = 0.55 ((c) 
and (d)). The period of the nodoid and height of the bubbles 
diverge when the two surfaces join. 
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4.1 Perturbative expansions 

Though an expression in terms of elliptic integrals is com- 
pact and useful for plotting, we have garnered further in- 
sight into the shapes via perturbation for small K. Start- 
ing with 



2z 
~R 



pudu 



(37) 



we expand to find this to second order in K: 
2z 




Away from r — 0, the first correction is negative, confirm- 
ing that the surfaces flatten out with increasing K. Note 
that the terms of order O(K) and above diverge at r = 0, 
indicating that the expansion is not uniformly convergent 
for all r. Examining (fl2)l . we see that the K — > limit is 
singular at r = since the term with the highest power of 
1/r drops out. 

The K — > oo limit does not suffer from this conver- 
gence issue. Expanding in powers of n = K , we find 



2z 
~R 



udu 



+ K 



u(l — u)du 

(1 - M 3 ) 3 / 2 



1 u 2 (l-u) 2 (-2-2u + u 2 )du 



2(1 - U 3)5/2 



(39) 



This expansion is regular for all < r < R for the first 
few terms, and one can argue that this perturbation series 
is regular in n and r since the singular point at K = 
0,r = has now moved to K — oo,r = 0. Fortunately, 
these integrals are related to known functions (useful for 
plotting!): 



2z _ VttT(|) 

R r(l) 



3 



/r\6 2 V 
\RJ ; 3'2, 




where B(z; a, b) denotes the incomplete Beta function and 
2^1 (a, 6; c; z) denotes the standard hypergeometric series. 
The first term gives us the K — > oo bubble shape, in other 



words the shape of a bubble with only nematic elastic 
energy and no surface tension energy. 

We have used these expansions to study the shape 
near the north and south poles at r = and the equa- 
tor at r = R. When K = the generating curve is a 
circle, quadratic near the poles and with a square-root 
cusp at the equator. However, for K ^ 0, the behavior 
is quartic at r = 0: z(r) — f^jra ?' 4 — 0(r 6 ). The power 
of K in the denominator in the coefficient of the leading 
term explains why the expansion in powers of K is sin- 
gular at r = 0. Near r = R, z(r) is still a square root: 

7 (f>_f\ _ (K+R 2 )V2R 

pansions are valid there. 



y/i-0{£ 3 / 2 ). Thus, both ex- 



5 Conclusion 

In the future, we will consider the dynamics of the two 
bubbles shown in Fig. [TJ The two limits we studied here 
perturbatively control the two long-time, scale-invariant 
limits of the solutions. As discussed earlier, at finite K, the 
problem has two length scales, R and \f~K - as the value of 

K = {^jf^j changes, one or the other limit becomes more 

appropriate, thus as the length scale set by geometry R 
(e.g. by assuming an initial volume) becomes much smaller 
than y/K through diffusion, we go to the K = oo or "pure 
nematic" limit. 

We were able to analyze this special case where the 
director was azimuthal and found the complete range of 
allowed shapes and topologies. Whether we can extend our 
method of analysis for longitudinal alignment or general 
alignment is an open question. 
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